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Résumé 

. On considère une classe de systèmes non-linéaires pour lesquels la 
technique de synthèse par backstepping n’est pas applicable. On présente 
un critère pour la synthèse d’une commande hybride faisant l’union d’une 
commande par backstepping et d’une commande locale. Cette commande 
hybride rend le système bouclé globalement asymptotiquement stable. Le 
critère de sélection repose sur le choix de la taille et la stabilisation d’un 
ensemble inclus dans le bassin d’attraction d’une commande locale. Les 
résultats sont illustrés par des simulations. 


1 Introduction 

Il existe de nombreuses méthodes différentes pour la synthèse de commandes 
pour les systèmes non-linéaires (voir par exemple, (Freeman and Kokotovié, 
2008; Khalil, 1992; Kokotovié, 1992)). En fonction de la structure de chaque sys¬ 
tème, la synthèse suit une approche différente, par exemple, avec un grand gain 
(Grognard et al., 1999), par forwarding (Jankovié et al., 1996; Mazenc and Praly, 
1996; Sepulchre et al., 1997) et par backstepping (Freeman and Kokotovié, 2008; 
Krstié et al., 1995; Praly et al., 1991). En raison de la présence des paramètres 
inconnus ou de la structure dynamique, certaines techniques peuvent être inap¬ 
plicables. 

Pour le type de systèmes où la technique de synthèse par backstepping peut 
ne pas être utilisable, l’approche présentée dans (Stein Shiromoto et al., 2011) 
donne une direction pour obtenir une commande stabilisante globale. Plus pré¬ 
cisément, la technique se propose de rassembler deux types de commandes. 
L’union de commandes ainsi que de fonctions de Lyapunov a déjà étudiée dans 
la littérature, voir (Andrieu et al., 2011), (Andrieu and Prieur, 2010) et (Prieur, 
2001). Dans cette dernière référence cette union a été réalisée en utilisant des 
commandes hybrides ((Hespanha et al., 1999), (Morin and Samson, 2000)). Un 

* E-mail: humberto. shiromotoOieee. org. 


1 



des aspects essentiels est que cette classe de commandes hybrides est robuste 
par rapport aux erreurs de mesure et aux erreurs d’actionneur. 

Il est important de noter que dans (Stein Shiromoto et al., 2011), une classe 
de systèmes non-linéaires a été présentée pour illustrer l’utilisation de cette 
technique. En revanche, cette approche laisse en suspens des questions ouvertes. 
Ainsi, il n’a été défini aucune méthode pour la détermination des gains des 
commandes obtenues. Dans cet article en reprenant l’exemple particulier étudié 
dans (Stein Shiromoto et al., 2011), nous montrons comment un choix approprié 
des paramètres nous permet de concevoir une loi de commande hybride qui 
stabilise globalement et asymptotiquement l’origine. Par ailleurs, on présente 
aussi des limitations sur cette technique, plus précisément on analyse un cas tel 
que l'hypothèse d’inclusion (voir hypothèse 3, introduite ultérieurement) n’est 
pas satisfaite. 

Dans la section 2 les concepts sur lesquels nous travaillons et la formulation 
du problème en considération seront introduits. Dans la section 3 nous présen¬ 
tons le résultat de ce travail pour la classe de systèmes non-linéaires considérée. 
Nous illustrons le résultat avec des simulations présentées dans la section 4. 

Notations. Dans cet article nous utilisons les notations suivantes : x ■ y est 
le produit interne cartésien entre deux vecteurs a: et y, || ■ || dénote la norme 
induite. La boule fermée et unitaire est notée par B. La classe des fonctions 
/ : R>o —> R>o qui sont continues, zéro à zéro, strictement croissantes et non 
bornées est notée par IC 00 . En ce qui concerne les dérivées, on note par d x f(x) 
la dérivée d’une fonction / par rapport à un vecteur x. La dérivée de Lie d’une 
fonction V par rapport au vecteur /, c’est-à-dire d x V{x) ■ f(x,u ) calculée en 
(x, u ) sera notée par LfV (x, u). Soit x = f(x, u(x)) un système en boucle fermée, 
son bassin d’attraction sera noté par B f(u). 

2 Formulation du problème 

Considérons le système non-linéaire défini par 

( Xi = fi{xi,x 2 ) + hi(xi,X 2 ,u) . . 

\ X 2 = f 2 (xi,x 2 )u + h 2 (xi,x 2 ,u) ' ‘ ' 

avec (. xi,x 2 ) G R" -1 x 1 et u € K la commande. Les fonctions /i, f 2 , h\ 
et h 2 sont localement Lipschitziennes. Par ailleurs, les fonctions décrivant la 
dynamique satisfont à l’origine /i(0,0) = h\ (0,0,0) = 0, h 2 (0,0,0) = 0 et 
f 2 {xi,x 2 ) £ 0, V(æi, æ 2 ) e R". 

On utilise aussi une notation plus compacte pour le système (1) : x = 
fh(x, u). Dans le cas hi = 0 et h 2 = 0, on note simplement x = f(x, u). 

2.1 Des hypothèses 

Hypothèse 1 [Stabilité locale] Il existe une fonction C 1 positive, définie et propre 
Ve, : R" —> R>o, une fonction continue tpg. : R™ —X R et une constante positive ce 
telles que, Væ G fl C( (Ve), 

Lf h V e (x,tpe) < 0, (2) 

avec Q ce (Ve) = {x : Ve(x) < c<, i / 0}. 
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Pour vérifier cette hypothèse, il est possible de faire la synthèse d’une com¬ 
mande locale pour (1) en exploitant une approximation du modèle. Par exemple, 
si l’approximation de premier ordre autour de l’origine est contrôlable, alors il 
est possible d’obtenir une commande linéaire locale et une fonction de Lyapunov 
quadratique telles que l’hypothèse 1 soit satisfaite. 

À cause de la présence des fonctions h\ et h 2 et leurs dépendances par rapport 
à la commande u, il n’est pas possible d’utiliser la technique de backstepping dans 
le sens classique. Plus précisément, supposons que l’on ait une commande ipi 
stabilisante pour le sous-système x\ = f(x 1 , £ 2 ) et une fonction de Lyapunov V\ 
associée (autrement dit que l'item a) de l'hypothèse 2 dessous est valide), alors, 
avec la procédure classique du backstepping la fonction V(x\,x 2 ) = V\(x\) + 
(x 2 — p\{x\)) 2 /2 est une fonction de Lyapunov candidate pour x = fh(x, P\(x))- 
On calcule V(x) := Lf h V(xi,x 2 ,u), V(xi,X 2 ,u) G R" -1 x R” x R, 


V(x) < -a(Vi(a;i)) + ( x 2 - ipi(xi)) • l/ 2 (xi,x 2 ) 
■u + h 2 {x 1 ,x 2 ,u) - Lf 1+hl <pi(x 1 ,x 2 ,u) 


+d Xx Vi(xi) 


/ fi(xi,sxi — (1 
J 0 


- s)^i(æi))ds] 


+L hl V 1 (xi,x 2 ,u). 


Soit 


E(xi,x 2 ,u) = [f 2 {x 1 ,x 2 )u + h 2 (x 1 ,x 2 ,u) 
-L fl+hl ip 1 (xi,x 2 ,u) 

+d Xl Vi(x 1 ) ■ / fi(xi, sxi — (1 — s)(/?i(æi)) ds]. 
J 0 


Alors, 


V(x) < -aiy^xi)) + (x 2 - ^i(a:i)) ^ 

■E(xi,x 2 ,u) + L hl Vi(xi,x 2 ,u) 

Par conséquence, pour obtenir un terme ( x 2 — ipi{x\)) 2 à droite de l’inégalité 
(3), il est nécessaire de résoudre le système d’équations implicites en u donné 
par 

L hl V 1 (xi,x 2 ,u) = -k(x 2 - ipi{xi)) 2 

E(xi,x 2 ,u) = -k(x 2 - tpi{xi)) 

avec k > 0. Ce système d’équations implicites est en général dificile à résoudre. 
Cette obstruction nous incite à introduire des hypothèses sur les fonctions h± et 
h 2 dans l’objectif de résoudre le problème de stabilisation. 

Hypothèse 2 [Bornes] Il existe une fonction V\ : R™ -1 —> R>o positive, définie, 
propre et de classe C 1 . Une fonction ip± : R" -1 —>• R de classe C 1 . Une fonction 
a : R>o — > R>o localement Lipschitz de classe K°°. Une fonction continue 
'L : R” —> R et deux constantes positives e < 1 et M telles que les propriétés 
ci-dessous sont valides 

a) (Contrôleur stabilisant pour le sous-système en x\) Vaq G R” - ' 1 , 

Lf 1 Vi(xi,(pi(xi)) < -a(Vi(xi)) 
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b) (Borne sur h\) M(x\,x 2 ,u) GW 1 1 x K x K, 

\\h 1 (xi,x 2 ,u)\\ < 'S>{x 1 ,x 2 ), 

L hl Vi(xi,<pi(xi),u) < (1 - e)a(Vi(x 1 )) + ea(M); 

c) (Borne sur d X2 h±) \/(xi,x 2 ,u) G R" -1 xRxI, 

\\d X2 hi(xi,x 2l u)\\ < ^(x\,x 2 )\ 

d) (Borne sur h 2 ) 'd(x\,x 2 , u) £ K" -1 x K x R, 

\\h 2 (xi,x 2 ,u)\\ < ^(xi,x 2 ) 

La condition a) de l’hypothèse 2 permet de concevoir une commande stabi¬ 
lisante globale pour l’équation définie par 


xi = fi(xi,x 2 ), (4) 

lorsque x 2 est vu comme une entrée. C’est-à-dire que cette condition permet 
de concevoir une commande pour x = f(x,u) à travers la technique de backs- 
tepping. Les conditions b)-d) fournissent des bornes pour les termes h\ et h 2 
ne permettant pas que cette technique soit appliquée pour le système (1). En 
exploitant les bornes sur les fonctions h\ et h 2l nous arrivons à concevoir une 
classe de commandes qui garantit que l’ensemble compact donné par 

A = {(xi,x 2 ) £ R" : Vi(xi) < M,x 2 = ip\(x\)}. (5) 


est globalement pratiquement asymptotiquement stable. Plus précisément, nous 
avons le résultat suivant présenté dans (Stein Shiromoto et al., 2011, Proposi¬ 
tion 3.1) 

Proposition 1 Si l’hypothèse 2 est satisfaite, alors l’ensemble A défini par (5) est 
globalement pratiquement stable, c’est-à-dire, pour chaque a > 0, il existe une 
commande continue <p g telle que A+aB contient un ensemble qui est globalement 
asymptotiquement stable pour le système x = fh{x, <p g (x)). 

La démonstration de la Proposition 1 fournit une commande ip g dont la 
définition est : 


V ‘ (n ’^= KvkL^' lKv 

L fl ipi(xi,x 2 ) - d Xl Vi(xi) 

■ / d X2 fi(xi,g Xl , X2 (s))ds 
J o 

~{x 2 - <pi{xi)) • (c + -A 2 (æi, æ 2 ))], 


( 6 ) 


avec Ky = 2 (M+a)/a 2 , r] XltX2 (s) = sx 2 + {l — s)ipi(xi), c une constante positive 
et 


A(xi,æ 2 ) 


||9 Xi El(xi)|| [ 'b(x 1 ,r] XuX2 (s))ds 
J o 

-\-'St(x 1 ,x 2 )K v (i + ||ô» æi <^i(m)||) 
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Afin d’obtenir une loi de commande qui stabilise globalement et asympto¬ 
tiquement l'origine, il faut donc rassembler la commande locale ipi donnée par 
l’hypothèse 1 et la commande obtenue à la Proposition 1 en exploitant un for¬ 
malisme hybride. Pour cela nous utilisons l’hypothèse suivante. 

Hypothèse 3 [Inclusion] La fonction Ve, satisfait l ’inégalité suivante 

max Ve(xi,x 2 ) < ce- 

(xi,x 2 )eA 

Les hypothèses 1 et 3 imposent que l’ensemble A soit inclus dans le bassin 
d’attraction de (1) en boucle fermée avec ipg. 

Dans (Stein Shiromoto et al., 2011), nous avons présenté le problème de 
conception d’une commande pour (1) de telle sorte que l’origine soit globa¬ 
lement asymptotiquement stable pour le système (1) en boucle fermée. Nous 
présentons ici la définition de ce type de commande (pour une lecture plus dé¬ 
taillée sur les systèmes hybrides, le lecteur est invité à consulter (Goebel et al., 
2009) et (Prieur et al, 2007)) 

Définition 2 Une commande hybride pour (1), notée par K, est composée de 

- Un ensemble totalement ordonné Q ; 

- Pour chaque q G Q 

- des ensembles fermés C q C R™ et D q C R" tels que C q U D q = R" ; 

- des fonctions continues ip q : C q —> R ; 

- des fonctions extérieurement semi-continues, localement bornées, uni¬ 
formément dans Q, évaluées en ensembles G q : D q =4 Q avec des images 
non vides. 


Avec ces données, nous pouvons définir une commande hybride telle que l’origine 
du système bouclé est globalement asymptotiquement stable (Stein Shiromoto et al., 
2011, théorème 1) : 


Théorème 1 Soient C( et ce deux constantes positives qui satisfont 0 < ce < ce ■ 
Si les hypothèses 1-3 sont toutes valides et s’il existe une constante a > 0 telle 
que la commande hybride K définie par Q = {1,2}, les sous-ensembles 


les commandes 


Cl = 

{{x\,x 2 ) G R" 1 x R : Vi(x 1 ,x 2 ) < q} 

c 2 = 

{{xi,x 2 ) G R" -1 x R : Vt(xi,x 2 ) > c/} 

D q = 

(I "- 1 

X R) \ Cq,Wq G Q, 

ipi : 

Ci 

->• R" -1 x R 


(xi,x 2 ) 

h> ipi(xi,x 2 ) = <fie(xi,x 2 ), 

■ 

c 2 

-» R” -1 x R 


(xi,x 2 ) 

ip 2 (Xi,X 2 ) = tp g (xi,X 2 ) 


et les applications multi-valuées et définies par D q 3 ( Xi,x 2 ) >—>■ G q (x 1 , 2 : 2 ) = 
{3 — q}, q G Q rend l’origine globalement asymptotiquement stable pour le sys¬ 
tème (1) en boucle fermée avec K ; 

( X = f h (x,<p q (x)), X G Cq 

\ q + g G q (x), x g D q 
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Remarque 1 Si les items b)-d) de l’hypothèse 2 sont valides, l’équation (3) té¬ 
moigne de l’effet des fonctions h\ et h^ sur la procédure de backstepping clas¬ 
sique. Intuitivement, plus les bornes des fonctions h\ et /12 sont grandes plus le 
bassin d’attraction de la commande locale est petit. Ainsi, il existe un compromis 
entre la taille des bornes des fonction hi, i = 1,2 et la validité de l’hypothèse 3. 

Une condition nécessaire pour que la commande hybride précédemment défi¬ 
nie rende l’origine globalement asymptotiquement stable pour le système bouclé 
est donnée par 

A Ç B fh (ip e ). 

Supposons que les hypothèses 1 et 2 sont valides. Si A % B f h (<pe)> alors l’hypo¬ 
thèse 3 n’est, pas satisfaite. Ainsi, il peut exister des conditions initiales telles que 
la solution du système x = fh{x,ip g (x)) converge vers l’ensemble A \ B f h (<pe). 

Ainsi les trajectoires du système x = fh{x,ipe(x)) ne convergera pas vers l’ori¬ 
gine pour toutes les conditions initiales dans l’ensemble A \ B f h (<pe). o 

Bien que le théorème 1 donne une commande globalement stabilisante, lorsque 
les hypothèses 1-3 sont valides, il n’est pas toujours évident d’effectuer la syn¬ 
thèse d’une commande locale qui satisfasse la contrainte sur le bassin d’attrac¬ 
tion. Par ailleurs, la taille de l’ensemble A dépend du gain de la commande 
stablisante pour le sous-système en x\. Dans ce travail, nous présentons des 
conditions sur les gains telles que l’hypothèse 3 soit satisfaite, quand la com¬ 
mande locale a été déjà donnée par l’hypothèse 1. 

2.2 Définition du problème 

Notre approche se décompose en deux étapes : 

Dans un premier temps, nous devons chercher un couple stabilisant, c’est- 
à-dire, une commande locale ipg et une fonction de Lyapunov locale Vt avec 
une estimation du bassin d’attraction pour le système x = fh(x, <pt(x)). 

Dans un second temps, nous devons établir des conditions sur le gain de 
la commande pour le sous-système défini par (4), telles que la taille de 
l’ensemble défini par (5) soit la plus petite possible afin que l’hypothèse 3 
soit plus facile à satisfaire. 

Dans cet article, nous n’allons pas analyser un problème général mais nous 
allons poursuivre l'étude d’une classe particulière de systèmes non-linéaires in¬ 
troduite dans (Stein Shiromoto et al., 2011). 

3 Une classe de systèmes non-linéaires 

On rappelle la classe de systèmes non-linéaires présentés dans (Stein Shiromoto et al., 
2011, Équation (12)) : 

/ *2 = u , 

\ x\ = X\ + Xi + 9x\ + 9(1 + X\) sin(u) v 

avec 9 > 0 qui est un paramètre libre mais constant pour chaque système. 

L’intérêt dans cette classe de systèmes nous permet d'illustrer l’influence 
de la norme des fonctions hi, i = 1,2, dans les calculs numériques et dans la 
faisabilité des hypothèses 1-3. 
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Dans (Stein Shiromoto et al., 2011, Lennna 4.1) nous avons vérifié que les 
hypothèses 1-3 sont valides pour tous les systèmes avec 0 < 0.001. Dans notre 
étude, nous montrons qu’avec un choix approprié de la commande <p\ les mêmes 
approches sont valides pour une classe plus grande de systèmes (plus précisément 
pour des valeurs plus grands de 0). 

L’idée de notre approche est de considérer un gain I\\ qui garantit que 
l’orientation du plus grand axe de l’ellipsoïde associé aux lignes de niveau de la 
fonction de Lyapunov soit celle de l’ensemble A. 

En suivant cette démarche, nous obtenons alors le résultat suivant. 

Théorème 2 Pour la classe de systèmes définie par 6 < 0.06, il existe une com¬ 
mande hybride telle que l'origine est globalement asymptotiquement stable pour 
chaque système en boucle fermée. 


Démonstration. 

Pour ne pas alourdir la lecture, les calculs intermédiaires seront omis, et nous 
ne présenterons que les éléments principaux de la démonstration. L’hypothèse 
1 est vérifiée, avec la commande 


<Pt{xi,x 2 ) = —X\k\ + x 2 k 2 


( 8 ) 


avec k\ = 7 + 0 et k 2 = —•4 + 40 + 0(1 + 9). La fonction de Lyapunov pour le 
système (7) en boucle fermée avec ipe est définie par 


Vé(xi,x 2 ) 


W T P(W = w T 


'5/2 

1 



vj = [xi — 0 x 2 ,x 2 ] t . La constante et vaut 


Q 


( 2 - 0 Pl { 0) \ 2 

\ 0p 2 {0) J ’ 


V0 < 01. 1 


(9) 


( 10 ) 


L’hypothèse 2 est vérifiée avec <fi(xi) = — 2.7456aq — 0x\, Vi(xi) = x\/2, 
a(s ) = 3.4912s, ^( xi,x 2 ) = 0(1 + |xi|), £ = 0.89 et M = 0.02. Par ailleurs, la 
proposition 1 assure l’existence de la commande donnée par (6), avec c = a = 10, 


tp g (x i,x 2 ) = 
x\ 

2 K v 


-(2.7456 + 20x 1 )(x 1 + x 2 + ôx\) 

X 2 -<Pl(X!) C 2 

--( C +4 A ( æ l> æ 2)), 


(H) 


et 


A(xux 2 ) = 6K V ( 1 + \xi\)(\xi\/K v + 1 

+ |2.7456+26tci|) 

telle que l’ensemble A = {(xi,x 2 ) : |xi| < 0.2, x 2 = ^i(xi)} est globalement 
pratiquement stable pour le système (7) en boucle fermée avec <p g . 

1. Les polynômes pi et p 2 sont donnés par 

pi (6>) = -396 + 23086» + 97686» 2 + 144O0 3 

p 2 (0) = 792 + 70006» + 208566» 2 + 216726» 3 + 21606» 4 . 

et d 1 = 0,118462 est la plus petite racine positive de l’équation £(6») = —2 + (ij> \ (0). 
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Pour l’hypothèse 3, on trouve qu’elle est vérifiée, Vd < 0.0607418. En parti¬ 
culier, avec 9 = 0.06 on a 

max Vt{xi,x 2 ) = Vé(0.2, y>i(0.2)) 

(xi ,æ 2 )£A 

= 0.030807 < 0.0390824 = c e 

Avec le théorème 1, nous obtenons donc une loi de commande hybride nous 
permettant de conclure la preuve du théorème 2. • 


4 Simulations 

Dans cette section, nous présentons des simulations du système précédant. 
Nous considérons comme condition initiale des points au bord d’une boule de 
rayon 2 et nous exécutons une simulation du système (7) avec 9 = 0.06 et bouclé 
avec K. 

Comme attendu, le système converge vers l’origine. Pour avoir une illustra¬ 
tion claire nous montrons l’évolution temporelle de Vu (figure 1) et des compo¬ 
santes de la solution (x, q) du système bouclé avec K seulement pour la condition 
initiale (x\,x 2 ,q) = (2,0,1) (figure 2). Pour les autres points qui appartiennent 
à la boule le comportement est similaire, comme le montre la figure 3 qui est la 
trajectoire de 5 points dans l’espace d’états. 



Figure 1 - Evolution de (9) pour le système (7) avec 9 = 0.06, en boucle fermée avec K. 
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Figure 2 - Evolution des composantes de la solution du système (7) avec 9 = 0.06, en boucle 
fermée avec K. 
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Figure 3 - Trajectoires des solutions pour différentes conditions initiales choisies au bord d'une 
boule de rayon 2. 
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5 Conclusion et perspectives 

Dans ce document, nous montrons que par un choix approprié de la loi de 
commande locale ipi pour le sous-système X\, cela permet de vérifier l’hypothèse 
3 pour une classe plus grande de systèmes que ce que nous avions obtenu lors de 
nos travaux précédents dans (Stein Shiromoto et al., 2011). A présent, l’objectif 
est donc de trouver une démarche constructive de sélection de ce contrôleur local 
afin de garantir que le bassin d’attraction de celui-ci contiendra l’ensemble A et 
par conséquent de vérifier l'hypothèse 3. 
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